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The proton-air inelastic cross section value o-^^^J^^^-^^=338±21{stat)±19{syst)-28{syst) mb at -y^ 2 TeV has 
been measured by the EAS-TOP Extensive Air Shower experiment. The absorption length of cosmic ray 
proton primaries cascades reaching the maximum development at the observation level is obtained from the flux 
attenuation for different zenith angles (i.e. atmospheric depths). The analysis, including the effects of the heavier 
primaries contribution and systematic uncertainties, is described. The experimental result is compared with 
different high energy interaction models and the relationships with the pp (pp) total cross section measurements 
are discussed. 



1. Introduction 

This energy region -^s « 2 TeV is of particular rele- 
vance because of high energy physics and astrophysics 
issues. 

The pp total cross section, cr*°*, and Cp'^fair ^'i'® i"®- 
lated and can be inferred from each other by means 
of the Glauber theory. The whole procedure is model 
dependent, the results [l|41| differing of about 20% for 
-y/s values in the TeV energy range. Available acceler- 
ators' measurements at the highest energies, are them- 
selves affected by systematic uncertainties of difficult 
evaluation. The pp (pp) cross section measurements 
at energies of = 1.8 TeV [l-Q differ of about 10%, 
exceeding the statistical uncertainties of the measure- 
ments. It is therefore of primary interest to have ex- 
perimental measurements of cr^ll^^ and cr*°* at the 
same center of mass energy, i.e. around y/s ~ 2 TeV, 
at which collider data are available. 

The interpretation of Extensive Air Shower mea- 
surements relies on simulations that use hadronic in- 
teraction models based on theoretically guided extrap- 
olations of the accelerator data obtained at lower en- 
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ergies and restricted to limited kinematical regions. A 
'^p^ih- direct measurement and the comparison of ob- 
servables as obtained from measurements and model 
based simulations, in the same conditions, is therefore 
highly recommended in order to confirm the validity 
of the whole analysis procedure. 

Measurements of the p-a,ir inelastic cross section 
performed in EAS have been reported. Since air 
shower detectors cannot observe the depth of the first 
interaction of the primary particle, indirect methods 
have to be used. Two main techniques are used: the 
constant Ng-N^ cuts [10- 12] by means of particle ar- 
rays, and the study of the shower longitudinal profiles 
using fluorescence detectors [l^[3| at higher energies. 

Following the particle array technique, the primary 
energy is first selected by means of the muon number 
(N^). Proton induced showers at the same develop- 
ment stage are then selected by means of the shower 
size (Ng). The cross section of primary particles is ob- 
tained by studying the absorption in the atmosphere 
(Aobs) of such showers, through their angular distri- 
bution at the observation level. The rate of showers 
decreases exponentially with zenith angle 9 (i.e. at- 
mospheric depth of first interaction) as: 

f{6) = G(0)/(O)expha;o(sec0- l)/Aobs] (1) 

where xq is the vertical atmospheric depth of the de- 
tector, and G{9) the angular acceptance. 

With air fluorescence detectors, the absorption 
length Xobs is obtained fltting the atmospheric depth 
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of the maximum shower development stage {Xmax) 
distribution taiL 

The observed absorption length in both cases is af- 
fected by the fluctuations in the longitudinal devel- 
opment of the cascades and in the detector response. 
Such fluctuations can be studied through simulations, 
providing the conversion factor k between the ob- 
served absorption length and the interaction length 
of primary protons {k=\obs / ^int) ■ This factor is then 
used to convert the observed experimental absorption 
length A^^^ into the interaction one A^^J'. 

In this paper we will report on the measurement 
of the p-air inelastic cross section at primary energy 
Eok2-W^ eV {i.e. Eo = (1.5 2.5)-10i5 eV, Vs » 
2 TeV) with the EAS-TOP experiment. Primary en- 
ergies are below the steepening {knee) of the primary 
spectrum (i.e. Eq < 3 • 10^^ eV) above which the 
proton flux is strongly reduced. The constant Ng- 
analysis has been optimized [ll| selecting show- 
ers at the maximum development stage where fluc- 
tuations are lower and heavier primaries rejection is 
improved by the request of the highest Ne values at a 
given primary energy. 

The constant Ng-N^ method with the selection of 
cascades' maximum developement stage is equivalent 
to the study of X^ax distribution tail. As shown in 
Fig. [U the accessible part of the X„iax distribution 
tail depends from the vertical atmospheric depth (xq) 
of the detector. This is a limitation on the possibility 
to exploit this method at different atmospheric depths 
and on the maximum zenith angle (i.e. atmospheric 
depth) that can be considered in the analysis without 
running out of statistics. Systematic uncertainties of 
the measurement and the effect of possible contribu- 
tions of heavier primaries are discussed and evaluated. 



2. The experiment and the simulation 

The EAS-TOP array was located at Campo Imper- 
atore, Gran Sasso National Laboratory, 2005 m a.s.L, 
820 g/cm^ atmospheric depth. 

The e.m. detector was made of 35 modules 10 m^ 
each of plastic scintillators, 4 cm thick, distributed 
over an area of 10^ m^. The trigger is fully efficient 
for Ng > 10^, i.e. for primary proton energies Eq > ?>■ 
10^^ eV. The experimental resolutions for Ng > 2 • 10^ 
are: ON^/^e — 0.1; gx^ = cry^ ~ 5 m; Cs — 0.1. 
The arrival direction of the shower is measured from 
the times of flight among the modules with resolution 
ae =i 0.9°. 

The muon-hadron detector (MHD), located at one 
edge of the e.m. array, is used, for the present analysis, 
as a tracking module with 9 active planes. Each plane 
includes two layers of streamer tubes (12 m length, 
3x3 cm^ section) and is shielded by 13 cm of iron. 
The total area of the detector is 12 x 12 m^. A muon 
track is defined by the alignment of at least 6 fired 
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Figure 1: Depth of shower maximum {Xmax) distribution 
for proton showers simulated with QGSJET II in the 
selected energy range {Eo = (1.5 2.5) • 10^^ eV) . 
The shaded area represents the interval of atmospheric 
depths (i.e. 1.0 < secO < 1.2) considered in the analysis. 

wires in different streamer tube layers leading to an 
energy threshold of Ejj^ « 1 GeV. 

A detailed description of the performance of the 
e.m. detector and of the muon-hadron detector can 
be found in Ref. 

EAS simulations are performed with the CORSIKA 
program 



17| with QGSJET 11.03 18] and SIBYLL 2.1 



[19| high energy hadronic interaction models These 
models have been widely employed for simulating at- 
mospheric shower developments and have shown to 
provide consistent descriptions of different shower pa- 
rameters in the considered energy range. Hadrons 
with energies below 80 GeV are treated with the 
GHEISHA 2002 interaction model. 

Proton and Helium showers have been simulated 
with an energy threshold of 10^^ eV, spectral index 
7 = 2.7 (7ffe = 2.65) up to 2 • lO^^ eV. KASCADE- 
like spectra 9] (resulting from our own fits) have been 
afterward sampled [lli] . 

Simulated events have been analyzed using the same 
procedure followed for experimental data. 



3. The analysis 

The rate of showers of given primary energy (Eq.i < 
Eo < Eo,2) selected through their muon number 
(N^.i < < N^.2) and shower size Ng correspond- 
ing to maximum development (Ng^i < Ng < Ne,2) is 
measured. 

The physical quantities required for the analysis are 
obtained through simulations, based on QGSJET II 
and SIBYLL interaction models, as described below. 
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Events in the specific proton primary energy range 
{Eq = (1.5 -f- 2.5) • 10^^ eV) are selected by means of a 
matrix of minimum (N^,i) and maximum (N^,2) de- 
tected muon numbers for every possible combination 
of zenith angle and core distance from the muon de- 
tector. The selection table is obtained from simulated 
data for 5 m bins in core distance (50 m < r < 100 
m) and 0.025 sec 6* bins (1.0 < sec6l < 1.2) for zenith 
angle. The selection of proton initiated cascades near 
maximum development is based on simulated distri- 
butions of the shower sizes at maximum development 
-^max |.]-^g sgigQ^;g(j energy interval. Choosing the 
shower size interval LogN™°-^ ± aLogN'""'' (i-e. 6.01< 
Log Ne < 6.17 for both interaction models) provides 
the selection of about 65% of the events around the 
maximum of the N™"^ distribution. 
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Figure 2: Acceptance corrected number of events vs. sec9 
for the simulated (solid circle for QGSJET II and solid 
triangle for SIBYLL) and experimental data selected 
with the N^-Ne cuts calculated with the two interaction 
models (open circle for QGSJET II and open triangle for 
SIBYLL). The fits with expression providing the Xots 
values are also shown (continuous lines). 

The interaction length Af™ is obtained as the av- 
erage proton interaction depth in the selected energy 
range, it results to be Af™ = 60.3 ±0.1 g/cm^ for 
QGSJET II and Af™ = 59.4±0.1 g/cm^ for SIBYLL. 

The acceptance corrected numbers of selected 
events N'^^i vs. zenith angle are shown in Fig. [51 
The fit with expression ([T]) provides A^^™ = 68.5 ± 
1.4 g/cm2 for QGSJET II and A^^^ = 69.9±1.4 g/cm^ 
for SIBYLL. Therefore k = A^/A™™ = 1.14 ± 0.02 
for QGSJET II and k = 1.18 ± 0.02 for SIBYLL 



4. Results 

The same analysis procedure discussed for the sim- 
ulations is applied to the experimental data. The 
corresponding event numbers as a function of sec{0) 
are also shown in Fig. [21 together with their fits 
providing XllP = 80.2 ± 4.3 g/cm^ and A^^f = 
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Figure 3: Shower size distribution (Log(Ne)) for proton 
(continuous) and helium (dotted) showers simulated with 
QGSJET II in the selected energy range . The two 
vertical dashed lines delimitate the Log(Ne) interval 
considered in the analysis. 

84.7 ± 5.0 g/cm^ for QGSJET II and SIBYLL re- 
spectively. From the relation ^int~K>bs/^^ obtain 
A^J = Ap_air = 70.7 ±4.2 g/cm2 for QGSJET II and 
A- J = Ap_air = 71.8 ± 4.5 g/cm2 for SIBYLL. The 
p-air inelastic cross section is then obtained from the 
relation crp^J^jj. (mb) = 2.41 • 10^/Ap_air, and results 
to be cr^"_°^i^ = 341 ± 20 mb with QGSJET II and 

4'-iir = 336 ± 21 mb with SIBYLL analysis. 

The contribution of heavier nuclei has been evalu- 
ated by simulating helium primaries with QGSJET II, 
assuming the KASCADE spectrum and composition, 
which accounts for an helium flux about twice that 
of the protons in the energy range of interest but, as 
shown in Fig. [31 due to the high Ng values requested in 
our analysis, the helium conatmination is reduced to 
less than 25%. The overall simulated observed absorp- 
tion length becomes \'^J-^(p+^'^'> — g2.6 ± 1.0 g/cm^, 
which implies k(P+^'=) = 1.04±0.02 , and cr^'!.!,^ 312 
±17 mb. Heavier primaries (i.e. CNO) hardly pass 
the N^-Ne cuts. 

The analysis procedure based on one interac- 
tion model has been applied to a simulated exper- 
imental data set produced with a different inter- 
action model with known p-air inelastic cross sec- 
tion in order to evaluate the systematic uncertain- 
ties. The result of the data simulated with SIBYLL 
(cr*'!f^i^=406±l mb) when analyzed with QGSJET 
II is cr;™^^;^=393±ll mb (Acr^™=^i,.=-13±ll mb) 
and viceversa CTp"'^^ij.=419±12 mb when the data 
simulated with QGSJET II (cT;'!.'=ii,=400±l mb 
(ACT™'=ii,=+19±12 mb) are analyzed with SIBYLL. 

An experimental data set has been simulated us- 
ing QGSJET II with HDPM ^ cross section value 
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(CTp'!f^;^=367±l mb) in order to check the capabil- 
ity to discriminate between two different values of 
the p-air cross section within the same interaction 
model. The lower cross section value is clearly dis- 
criminated by the anlysis based on SIBYLL that gives 



-5±13 mb). The dif- 



a»l,=372±13 mb (Aa»l,- 
ferences between the simulated and measured values 
(AtTpl'^j^jj.) are both positive and negative and compat- 
ible with the statistical uncertainties. Therefore we 
define as maximum systematic uncertainty the value 
(^syst = 19 mb which provides the value corre- 
sponding to 99% C.L for the distribution of the four 
deviations Aap"fl^j.. 
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Figure 4: p-air inelastic cross section data including the 
present measurement (solid square), and different 
hadronic interaction models. 



5. Conclusions 

Combining the results obtained with the two consid- 
ered interaction models and including the systematic 
uncertainties the |?-air inelastic cross section is: 



p— air 



= 338 ± 2\,,„,t ± 19 



syst 



As shown in Fig. 21 this value is about 15% smaller 
than the values in use within QGSJET II and SIBYLL 
and in better agreement with Refs. 0, [13, |2l|. Pre- 



dieted 



values, obtained from different (7^°' Teva- 



tron measurements at ^/s — 1.8 TeV by using dif- 
ferent calculations based on the Glauber theory, are 
reported in Fig. [SJ The present measurement is con- 
sistent with smaller values of the pp total cross section 
(ct*;*=72.8±3.1 mb [tJ, and ct*;*=71±2 mb |]), and 
the pp to p-air calculations predicting for a given value 
of CTpp* , a smaller value of af^ '^^ H, . Independently 
from the cross section analysis, the measured values 
of the absorption length Xots are about 15% higher 
than the simulated ones for both the considered in- 
teraction models. This can be probably ascribed to 
the fact that the measured cascades penetrate deeper 
into the atmosphere than predicted by the interaction 
models. 
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Figure 5: p-air inelastic vs. pp total cross section data. 
The present result (± 1 s.d., solid lines) is shown 
together with the results of different calculationsderived 
from pp measurements reported at ^/s = 1.8 TeV. 
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